arXiv:1506.08057v2 [cond-mat.quant-gas] 12 Sep 2015 


Methods for detecting charge fractionalization and 
winding numbers in an interacting fermionic ladder 


Leonardo Mazza\ Monika Aidelsburger^’^, Hong-Hao Tu^, 
Nathan Goldman^’^, Michele Burrello^ 

^NEST, Scuola Normale Superiore and Istituto Nanoscienze-CNR, Piazza dei 
Cavalieri 7, 1-56126 Pisa, Italy 

^Max-Planck-Institut fiir Quantenoptik, Hans-Kopfermann-Str. 1, D-85748 Garching, 
Germany 

^Fakultat fiir Physik, Ludwig-Maximilians-Universitat, Schellingstr. 4, D-80799 
Miinchen, Germany 

■^Center for Nonlinear Phenomena and Gomplex Systems, Universite Libre de 
Bruxelles, GP 231, Campus Plaine, B-1050 Brussels, Belgium 
^Laboratoire Kastler Brossel, College de France, 11 place Marcelin Berthelot, 
F-75005, Paris, France 

Abstract. 

We consider a spin-1/2 fermionic ladder with spin-orbit coupling and a 
perpendicular magnetic field, which shares important similarities with topological 
superconducting wires. We fully characterize the symmetry-protected topological 
phase of this ladder through the identification of fractionalized edge modes and non¬ 
trivial spin winding numbers. We propose an experimental scheme to engineer such a 
ladder system with cold atoms in optical lattices, and we present two protocols that can 
be used to extract the topological signatures from density and momentum-distribution 
measurements. We then consider the presence of interactions and discuss the effects 
of a contact on-site repulsion on the topological phase. We find that such interactions 
could enhance the extension of the topological phase in certain parameters regimes. 
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1. Introduction 

The experimental engineering of topological phases of matter in ultracold atomic 
gases [D El El la El E] lays the foundations for a deeper understanding of phase 
transitions that transcend the Landau paradigm of symmetry breaking. In these 
experiments, models displaying non-local forms of order are realized in highly 
controllable environments, where the parameters driving the system in and out 
the topological phases can be tuned with wide freedom, and where observables 
complementary to those of a typical solid-state experiment can be measured. These 
experiments have already revealed interesting properties associated with 2D topological 
Bloch bands: the anomalous (Hall-like) velocity, which was detected in response to 
an external force Hi], the topologically invariant Chern number [1], and chiral edge 
currents mum- 

Following these advances, an important objective would be to probe the edge 
modes of ID topological systems, which typically appear at zero energy and exhibit 
charge fractionalization. In particular, identifying an observable that unambiguously 
signals their presence in experiments constitutes a remarkable challenge. Detecting 
the properties of zero-energy edge modes would strongly complement the Zak-phase 
measurement, recently demonstrated with bosonic atoms in a ID optical superlattice 

[I]- 

Several theoretical efforts have been devoted to the design of realistic platforms 
hosting topological superconducting phases with ultracold fermions m El El uni El in 
nE], including number-conserving setups mminiEiiiHiiisiEQ]. Here, we envision 
an even simpler scenario, based on the fact that similar topological edge physics can 
be accessed without pairing mechanisms. Indeed, it is a generally overlooked fact that 
pairing interactions are not strictly necessary to mimic topological superconductors. 
A fundamental example is offered by the Su-Schrieffer-Heeger model IZII. which 
presents particle-hole symmetry and belongs to a non-trivial topological class of chiral 
Hamiltonians, namely the class BDI of the Altland-Zirnbauer classihcation [221 ESI EH 
EH EE]. Similar one-dimensional fermionic models without superconducting interactions 
display topologically-protected edge modes localized at their boundary, which are Dirac- 
like [27] and feature remarkable properties, such as charge fractionalization [28] . 

The goal of this article is twofold. First, we propose a route to mimic the physics of 
topological superconductors using state-of-the-art ultracold fermionic experiments. We 
exploit a simple idea: employing a two-leg ladder to double the fermionic species, in such 
a way that the legs are respectively associated with effective holes and particles [221EO] , 
see Fig. This geometry is indeed well-suited for experiments on ultracold gases, 
as it has already been realized for atoms trapped in optical lattices, either in 
physical ladder geometries [2], or exploiting internal degrees of freedom as an artihcial 
dimension [3 El ED EH E3]. We characterize the topological properties of this model, 
both in the absence and presence of interactions, hnding several affinities with previous 
studies of ID topological interacting fermionic systems [33l EH EH [Ml ED EH EH HQ]- 
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We argue that such systems constitute a useful playground, not only to examine the 
appearance of symmetry-protected edge modes, but also to study the role of contact 
interactions, which may be tuned to drive transitions between trivial and topological 
phases. This analysis of the Hubbard repulsion extends to spin 1/2 fermions previous 
studies about the effect of interaction in generalizations of the SSH model HD Eg. 

As a second objective, we analyze in detail how topological signatures might be 
directly observed in such interacting systems. Beside the detection of fractionalized 
edge-modes, we also focus on the winding number associated with the expectation value 
of the spin, which provides a good detection tool for topological phases also in the 
presence of interactions. In this way, we extend to a one-dimensional model in the 
topological class BDI the techniques developed to reveal the topology of cold-atom 
realizations of the two-dimensional Haldane model El] and other two-dimensional 
topological systems ESI EHl EZl EH]. For both these observables, we examine the effect 
of a trapping potential, which sets soft boundaries to the system (usually considered to 
alter the observation of edge physics). Our schemes are based on the direct observation 
of the atomic cloud or on time-of-flight measurements: in both cases they are extremely 
robust to such confinement. 

This article is organized as follows. In Sec.[^we introduce the model and provide an 
intuitive description of its symmetries. In Sec. [^we focus on the non-interacting model 
at half-hlling and thoroughly characterize the topological insulator that is reached for a 
certain range of parameters. Furthermore, methods to detect unambiguous signatures of 
the topological properties are proposed, based on the density profile and the momentum- 
distribution of the gas. In Sec. we study the role of interactions and characterize 
the related interacting topological insulator. In Sec. we describe a possible physical 
realization of the model based on laser-assisted tunnelings and in Sec. we present 
our conclusions. Finally Appendix A presents a detailed analysis of the non-interacting 
model, its order parameter and spin winding number. 


2. The model 

We consider a spinful fermionic ladder in the presence of external gauge potentials as 
depicted in Fig. [Tj Here the two legs of the ladder are associated with a pseudo-spin 
and the lattice sites along the main axis (x direction) are labelled by r = 1, ...,L. We 
introduce the 4-component fermionic operator a^, defined on the lattice site r, which 
acts on both the pseudo-spin and the spin related to the two internal states of the 
fermion (two commuting sets of Pauli matrices and are used to describe these 
degrees of freedom). The Hamiltonian describing the ladder system is taken to be of 
the form {h = 1) 

Ho = ^ h.c. j -f dl + Jt^ + + fio) hr} , (1) 

r 

and is schematically represented in Fig. The first term in Eq. ([^ describes the intra¬ 
chain tunneling along the x direction, with hopping amplitude t and spin-dependent 
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Figure 1. Schematic representation of the non-interacting model. The system is 
engineered in such a way that the two chains in the ladder present opposite kinetic 
energies. This is obtained through the introduction of a tt flux in each ladder plaquette. 
Here, the pseudo-spin t refers to the two legs of the ladder, while the spin tr is associated 
with the internal states (“spin”) of the atoms, see Eq. Q. 


Peierls phase-factor exp{iBaz/2), which represents a “spin-orbit coupling” analogous to 
those already realized in fermionic j5] and bosonic [6] chains. Note that, because of the 
Tz factor, the two chains have opposite kinetic energy; consequently, the motion around 
each plaquette of the ladder acquires a vr-phase independently of the cTz-spin component. 
The second term describes an on-site spin-flip term with amplitude hi, the inter-chain 
tunneling with amplitude J, the potential difference between the two chains fi, and the 
overall chemical potential fio- 

In the following we will elaborate on the fact that this system effectively reproduces 
some of the physical features of topological superconductors (TSC) without recurring 
to any physical pairing mechanism. This analogy is based on the idea that atoms 
in the first chain (r^ = -f-l) can be identified with conduction electrons of a generic 
ID superconducting model, whereas those in the second chain 

{tz = —1) can be identified with its holes, = cjjoie(?")- In this picture any 

tunneling from one chain to the other constitutes an effective pairing interaction, i.e. 

= cljg^^Choie ~ ^lect^Lct) where the last equality is justihed by the 
Bogoliubov-De Gennes treatment of the superconductor. Within this parallelism we 
interpret the four-band Hamiltonian ([^ as a Bogoliubov-de Gennes Hamiltonian in the 
superconducting picture. Specihcally, the particle-hole symmetry, which plays a key 
role in the physics of TSGs, is here represented by a swap of the two chains, C = Tyay, 
which have opposite kinetic energy in the same way as particles and holes do. Such 
mapping, though, must be seen only as an analogy, since the number of degrees of 
freedoms in the system ([^ is doubled with respect to the superconducting wire and this 
has important physical consequences, as will be discussed in the following. Finally, note 
that Hamiltonian ([^ is unitarily related to those considered in the four-wire setup of 
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Refs. |29l ISnj and may have a relevance also for the stndy of electronic gases. 

The model in Eq. Q can be realized using cold atoms trapped in an optical 
lattice. We present here an overview of the experimental proposal and refer the 
interested reader to Sec. for a detailed analysis of the implementation of the 
model. Let’s start considering a two-dimensional setup. The realization of a spin- 
dependent intra-chain tunneling, described by the hrst term in Eq. Q, is particularly 
challenging, as it requires a subtle control over the hopping amplitudes. This effect 
could be engineered by exploiting the laser-induced-tunneling methods implemented 
in recent experiments l50] . Specihcally, we propose to achieve this task by 

combining a spin-dependent staggered potential with large energy offset A between 
neighboring sites, inhibiting the bare hopping along the x direction, together with an 
onsite energy modulation set at the resonant frequency u = A. The spin-dependent 
staggered potential is chosen to be opposite for the two internal states, i.e. 14tag(?") = 
(— 1 )'’(A/ 2 )(T 2 , which can be realized by considering an appropriate anti-magic wave¬ 
length this choice is motivated by the fact that the resonant modulation will then 
generate effective tunneling matrix elements of the desired form tesi^) = t exp[i(j){r)az], 
see Refs. [IHl 1S21 ES]- In order to make the Peierls phase-factors constant over the 
whole lattice, i.e. 0(r) = R/2, we propose to modulate the lattice with two pairs of 
lasers; such a conhguration allows to address individual links independently |1], hence 
realizing the desired Peierls phase factors on all links, (see Section]^. Using additional 
helds resonant with the energy difference between the two spin-states a tunable onsite 
spin-flip term VlCx can be realized. Finally, the two-leg ladder can be isolated using an 
additional superlattice, or a light-intensity mask [511E5]. 


3. Topological phases in the non-interacting system 


Hamiltonian ([^ is characterized by four energy bands; as shown in Fig. by varying 
the hlling of the ladder, and thus the chemical potential /xq, the system is driven through 
a series of metallic and insulating quantum phases (we consider in this article only the 
case of zero temperature). For half hlling, corresponding to the case where /xq = 0 
and the particle density is p = N/L = 2 {N is the number of fermions), the single¬ 
particle Hamiltonian shows both the particle-hole symmetry we sought for, dehned by 
the operator C = Tyay and an additional time-reversal symmetry, T = which bring 
the system into the topological class BDI (see Appendix A for more details). This 


class includes, for example, the Su-Schrieffer-Heeger (SSH) model and, according to 
the periodic table of topological insulators and superconductors [2H123] , h may present 
topological phases with zero-energy modes. Specihcally, our model displays a non-trivial 
topological insulating phase for < ^c, 2 i where are dehned, for B < n, as: 


= r + {p-2cos{B/2)y 


+ {p + 2cos{B/2)y 


( 2 ) 


We hnd that the topological phase is surrounded by two topologically trivial phases. For 
Po = 0, trivial and topological phases are distinguished by a topological order parameter 
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Figure 2. Phase diagram of Hamiltonian Q. (Left): Energy of the eigenmodes 
of Hamiltonian 0 for periodic boundary conditions (red lower triangles) and open 
boundary conditions (blue upper triangles) for the parameters B = 7r/2, J/t = 1.0, 
il/t = 1.75, g/t = 0.5, Ufj/t = 0; the system size is L = 200. The lower-right inset 
zooms into the zero-energy region and shows the existence of two zero-energy modes for 
the open system. The upper-left inset shows the squared modulus of the wavefunctions 
of these two modes, which are localised at the edges. (Right): The phase diagram as a 
function of g^/t is derived from the previous spectrum: it alternates between normal 
insulating phases (Nl) and metallic ones (M). For half filling, thus at density p = 2 
{go = 0), the system is in the topological regime (BDI). 


W which takes the respective values of +1 and —1 (see the Appendix and in particular 
Eq. (A.7) for the dehnition which is based on the technique developed in 


3.1. Fractionalized edge modes 

In the topological phase, two zero-energy fermionic modes appear in ladders with open 
boundary conditions, as showed in Fig. These modes are exponentially localized 
at the ends of the system and have important consequences on the density distribution 
characterizing the topological insulating phase when N = 2L + 1 fermions are introduced 
in the system with hard-wall boundary conditions, as displayed in Fig. We observe 
that such modes are described by Dirac operators and they are not Majorana modes 
as it would be expected in the superconducting analog wire. The hgure shows that 
the density in the bulk of the system indeed corresponds to the expected value p = 2. 
Moreover, analogously to the SSH model, a charge 1/2 is exponentially localized at 
each boundary. This important signature of charge fractionalization can be suitably 
identihed through the expectation value of the operator n* = J2m=i('^rn — 2) where 
hm = dj^dm, which detects the excess density with respect to the bulk value p = 2. As 
illustrated in the insets of Fig. an overall excess density of 1/2 is localised within a 
few sites from the left and right edges of the sample. 

Such a signature can be observed even in the presence of a harmonic conhnement. 
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Figure 3. Charge fractionalization and spin winding number in a system with hard- 
walls boundary conditions. (Left): Density profile of the system with N = 2L + \ 
particles. The simulation parameters are the same as in Fig. The inset displays the 
expectation value of n*. (Right): top and middle panels, (£(fc)) for a system with 2L 
fermions and with L fermions; bottom panel, S2{k)/2 for L fermions artificially loaded 
into the second band, see Eq. (§. The value W can be extracted only in the latter 
case. 


described by the following contribution to the Hamiltonian: Htr = 'Yfr'^rofdr, with 
Wr = w{r — L/2y. The effect of an external potential can be understood in a Thomas- 
Fermi approach as a space-dependent chemical potential fio{r). Due to the four energy 
bands, the system has three insulating phases for intermediate Filings, and, in the 
presence of a harmonic trap, this yields a typical wedding-cake structure with integer 
density plateaus (see Fig. hrst column). Remarkably, even in the presence of the 
trapping it is possible to identify the fractionalized modes, as we see next. These zero- 
energy modes extend in the intermediate metallic region between one trivial plateau 
(p = 1,3) and the topological one (p = 2), up to exponential corrections. 


In Fig. 


second column, we show the expectation value of n** = J2m=o('^L/2+m — ^), 
which is particularly suited for the detection of fractionalized edge modes in cases where 
the density in the center of the trap is p = 2. Moving from the center of the trap to 
the next plateau (p = 1), this operator measures the excess density with respect to the 
p = 1 value. One can obtain either an integer (no fractional modes) or an half-integer 
value (presence of one fractional mode). This is an unambiguous signature of the non¬ 
trivial topological phase (third row): indeed, in this case, the quantity {h**) becomes 
half-integer for values of j corresponding to the distance of the p = 1 plateau from the 
center of the trap. 

As a hnal remark, let us stress that the problem of detecting fractionalized edge 
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Figure 4. Charge fractionalization and winding number in a system with harmonic 
confinement for topological and non-topological phases. We consider a trap with 
w/t = 0.03 for B = tt/ 2, J/t = 1 and g/t = 0.5, for which ric,i/t — 1-35 and 
ilc, 2 /t — 2.16. Different rows refer to different values of fl/t, from up to down, 0.85, 
1.1, 1.75, 2.5. The number of particles is chosen to have a density p = 2 in the center of 
the trap. The first column displays the density of the system (hj). The second column 
displays (h**). The third column shows the winding number relative to the second 
band S 2 {k) computed for a system with trap (black) and without a trap (blue). 


modes through a density measurement was first addressed in Ref. ISZIEH], where this 
detection method relies on the optical measurement of reflected light. The recent 
experimental advances, however, allow for the challenging method presented above, since 
the feasibility of a single-atom detection for ultracold fermions in optical lattices has 
indeed been demonstrated [59l [60l |6T] . In particular, a combination of laser cooling and 
fluorescence detection enables an unambiguous measurement of the occupancy of single 
sites for both [59l [60] and ®Li [6l] gases. This is of particular importance for the 
detection scheme that we are proposing, because it could suffer from the experimental 
inability to fix the total number of atoms which are used in the many experimental 
realizations necessary to reconstruct the signal {nj). The novel single-atom microscopes 
will also allow a post-selection based on the global number of particles of the system, 
necessary to obtain an accurate measurement. 
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Another interesting signature of the topological phase is offered by the expectation value 
of the spin operator: 

t,{k) = ^dl{(T ®To)dk (3) 

where dk is the four-component annihilation operator in momentum space. This 
approach is inspired by the techniques presented for two-dimensional systems in 
Refs. [Ill HU uni 051 07] and for ladders in Ref. [62], where it was shown that 
the expectation value of the spin (S(A;)) provides a good observable to identify 
the topological invariant (winding number) of certain topological insulators. In the 
following, we generalize this procedure to our quasi-one-dimensional ladder model and 
show that also here the topological invariant W, which clearly identifies the non-trivial 
topological regime, can be extracted from (S(fc)). This is thus another example of the 
interesting concept that a time-of-flight measurement can detect topological order. 

The Hamiltonian ([^ can be written in a real form thanks to its symmetries; 
therefore {'Ey{k)) = 0 for each eigenstate of the system and (S(fc)) always lies in the 
X — z plane. The “spin winding number” S is defined as the number of times the vector 
(I](/c)) encircles the origin for k going from 0 to 27r (the lattice spacing is set to 1). Let 
us denote with S 2 {k )/2 the expectation value of the spin operator ([^ for a state of non¬ 
interacting fermions hlling completely and solely the second energy band. Remarkably, 
the parity of S 2 , which is the winding number of S 2 {k), coincides with the topological 
order parameter of the model: 

W = (-1)'^^ (4) 


see 


Appendix A for a demonstration and for details on the analytical calculation of 


this topological index for this specific model). Eq. (|^ is analogous to those derived for 
several other two-dimensional models [131011,05106107]: it relates a topological invariant 
to a quantity, S 2 , to be extracted via time-of-flight imaging. A similar behavior was 
discussed in Ref. [62] for a two-band generalizations of the SSH model. In our case, 
we stress that the second band of the model is the lowest-energy band with non-trivial 
topological order; this is also related to the fact that the fractionalized edge modes 
appear in the second bulk gap. 

We now describe how to measure S 2 , for realistic systems, even in the presence of 
a harmonic trap. The main problem is that the spin winding number has to be probed 
for the second band of the Hamiltonian only: in a physical realization of the topological 
phase, both the first and the second band are filled, and the acquired signal includes 
information of both. The right column of Fig. shows the expectation value of the spin 
(S(/c)) obtained when the system with hard-wall boundaries is hlled with 2L fermions 
(top panel) and with only L fermions (middle panel). Both signals are not particularly 
interesting. If we consider the artificial situation where atoms populate the second 
band only (bottom panel), the spin expectation value is characterized by a winding 
number that reproduces the behavior of W and encircles the origin in the topological 
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phase. In the ideal case of a hard-wall conhning potential, the required value of the 
second band can be extracted by repeating the experiment twice, at densities p = 1 
and p = 2: the difference of the measured distributions returns the sought information 
(S(fc))p=2 - (S(fc))p=i = 

In the presence of a harmonic trap, the wedding cake density profile suggests that 
the many-body wavefunction can be roughly thought as a state where each energy band 
a is uniformly populated by atoms (A^i > N 2 > > tf). In this case the 

measurement of the observable (S(A;)) returns: 

{tm = ( 5 ) 


where Sa{k) is the expectation value of the spin calculated in the thermodynamic 
limit for the single particle eigenstate of the energy band (see Appendix A). C 


is the discretization adopted for the Brillouin zone in the time-of-flight imaging (see, for 
example, ^31 H6] ). 

If we consider the case in which the density profile shows only two plateaus, the 
value of S 2 {k) can be estimated by comparing the observed (S(fc))p =2 with that of a 
realization with a single plateau only, {Y,{k))p=i\ 


S^{k) = ^ ( 6 ) 

where Ni and N 2 are the occupations of the two bands for the state with two plateaus, 
and A^( is the total number of atoms in the reference state with a single plateau. All 
the quantities A^i,A "2 and A^( can be experimentally accessed and we report in the 
right column of Fig. |^the comparison of the data obtained for hard wall and harmonic 
potentials. Our numerical simulations confirm that even in the presence of the trap S 2 
is equal to ±1 in the topological phase, whereas in the trivial phases, it is either 0 or 
±2 (see Fig. 1^. 

Let us conclude with some information on how to measure (S(A;)) through spin- 
resolved time-of-flight imaging |13l HU HHl SSI SZj in our setup. Special care is required in 
time-modulated systems with spin-dependent features |63], as considered in the specific 
proposal detailed in Sec. [^because spin-dependent observables can potentially undergo 
large and complicated micro-motion (rapid motion with a time-scale of the order of the 
driving period 2tt/u}), which typically alters the accuracy of measurements. In such 
schemes, stroboscopic measurements performed at specific times, 27r/ci; x n where n 
is integer, are generally required to extract relevant information relative to the spin- 
dependent quantities [HSI- For the scheme detailed in Sec. the micro-motion can be 
estimated from the unitary operators K{t) and R{t) dehned in that Section, through the 
method of Ref. [S3]. We find that (S^(fc)) is unaffected by the micro-motion; in contrast, 
an accurate analysis of (Sa;) does require a stroboscopic measurement. Moreover, we 
note that measuring the expectation value of Sj, also necessitates a 7r/2 pulse, which has 
to be short compared to the driving period in order to probe the system stroboscopically. 
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Figure 5. Phase diagram of the interacting model at half-filling in the U/t and 
J/t plane. (Left) Single-particle gap Ai/t. (right) Two-particle gap A 2 /t. The 
topological region is characterized by Ai = 0 (dark blue region in the left panel) 
and it is delimited by gapless regions defined by A 2 = 0 (dark blue regions in the 
right panel) as represented schematically in Fig. The calculations are performed for 
Vl/t = 1.8, B = f and g/t = 1 at L = 72 with bond dimension D = 200. 


4. Interacting system 

Let us now consider the role of interactions, with a special emphasis on the robustness 
of symmetry-protected topological order. It is experimentally relevant to consider an 
on-site Hubbard interaction in each leg: 

H'lat ^ 'y 1 0 -^= 10 - 2=1 • ('^) 

r,Tz 


and to analyze the phase diagram of Ho+Hi^t at half hlling, p = 2, which is characterized 
by the competition between the topological insulator (TI) and Mott insulator (MI) 
occurring in the presence of a strong contact repulsion. We employ a density-matrix 
renormalization group algorithm based on a Matrix-Product State (MPS) ansatz [6^165]. 
We will consider systems with open-boundary conditions with L = 72 and maximal bond 
dimension D = 200. 

The transition between TI and MI can be located via the analysis of the charge 
gap aX N = 2L. In particular, the single-particle gap is dehned as: 

/\i{N) = E{N+ l) + E{N-l)-2E{N), ( 8 ) 

where E{N) is the ground-state energy of the system with N fermions. Clearly, 
Ai(Af) > 0 for the MI because the system has a thermodynamic gap. On the other 
hand, the TI has zero-energy modes which ensure that E{N — 1) = E{N) = E{N + 1) 
and thus Ai(iV) = 0. Unfortunately, the mere calculation of Ai does not permit to 
discriminate the TI from a generic gapless phase, for which Ai(iV) = 0 too. We thus 
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Figure 6. Schematic structure of the phase diagram at half filling as a function of J 
and U. (Left) The blue region represents the topological region as extracted from the 
numerical results presented in Fig. (note that strictly speaking our numerics could 
not access the thermodynamic limit). (Right) Qualitative extrapolation of the phase 
diagram. Three gapped phases can be detected: topological insulator (TI), trivial band 
insulator (BI) and Mott insulator (MI). Blue lines represent the phase transitions. Our 
numerical results suggest that the MI and BI phases are adiabatically connected. The 
exact nature of the phase diagram in the dotted region cannot be established due to 
numerical limitations. 


consider also the two-particle gap: 


/\2{N) = E{N + 2) + E{N - 2) - 2E{N). 


(9) 


Whereas for a gapless phase A 2 (A^) is also equal to zero, for a TI it is larger than zero, 
signaling the gap which is protecting the phase. 

Based on this discussion, we now consider a systematic study of the Hamiltonian, 
focusing on the competition of the two terms which are responsible for a gap opening, 
namely, the interaction term proportional to U responsible for the MI, and the interchain 
tunneling proportional to J. Roughly speaking, we identify the pairing term as 
the one inducing the TI, since at f/ = 0 the system is in a topological phase for 
< 


0 


< = Q2 _ (2cos(R/2) - fxf 


see 


Appendix A for more detail) and the 


two chains decouple at J = 0. 

Fig. [^presents the numerical results for Ai(2L) and A2(2L) in the parameter space 
spanned hy U/t and J/t. The other parameters are chosen such that at (7 = 0 there is a 
TI, and are listed here for completeness: fl/t = 1.8, B = and fi/t = 1. Calculations 
are limited to the size L = 72 and a systematic scaling to the thermodynamic limit, 
as well as the exact evaluation of the properties of the critical lines, is beyond our 
numerical possibilities; additionally, the two-dimensional space is studied with a grid 
of 0.2 along both axis. Despite these limitations, the qualitative nature of the phase 
diagram emerges quite clearly. Indeed, through the study of Ai and A 2 we are able 
to identify the TI, the MI and the critical regions which separate them, resulting in 
the schematic phase diagram presented in Fig. The topological region is identihed 
with the large region where Ai = 0 but A 2 > 0 whereas for the MI both Ai and A 2 
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U/t t/U 


Figure 7. Energy differences of the ground states of the system for several fillings 
around N = 2L. The same data are plotted as a function of U/t (left) and t/U (right). 
The Hamiltonian parameters are J/t = 1, H/t = 1.8, B = f so that the cut 
corresponds to the line J/t = 1 in the phase diagram of Fig. L = 72 and the 
maximal MPS bond dimension is I? = 150 

are larger than zero. Critical regions with Ai = A 2 = 0 separate the two insulators. 
The schematic phase diagram discriminates the MI, whose appearance is driven by the 
on-site repulsion, from the trivial band insulator (BI), which appears instead also at 
t/ = 0. Our investigation did not identify a phase transition between these two trivial 
insulating phases, which are adiabatically connected. 

To better analyze the transition between TI and MI, in Fig. [^we focus on the line 
at J/f = 1, which entails a phase transition for U/t = Ucr ~ 1-9. We show the behaviour 
of the chemical potentials E{2L + a) — E{2L + a — \) for a = +2, -|-1,0, —1 as a function 
of U. Two qualitatively different behaviours are separated by Ucr- For U/t > Ucr 
the energy cost for adding one particle to the states with N = 2L or N = 2L + 1 
is approximately U (especially for large values of U/t). Conversely, subtracting one 
particle from the states with N = 2L or N = 2L — 1 does not yield any energy gain. 
Thus, Ai{N = 2L) and A 2 (A^ = 2L) are larger than zero and are approximately equal 
to U and 2[/, respectively: these are typical signatures of a MI. 

For U/t < Ucr the energy cost and gain for adding and removing one particle 
to/from the state with N = 2L are both equal to U/2. We interpret this as a signature 
of the fractionalization of the zero-energy modes of the TI: the charge excess n 1/2 
on top of the density plateau p = 2 (one particle per site) does cost a repulsive energy 
Un. Since the zero-energy modes have fermionic nature, they cannot accommodate 
more than one particle: the energy cost for adding one additional particle to the state 
with N = 2L + 1 becomes significantly larger than U/2 (vice versa for removing one 
particle from the state with N = 2L — 1). Thus, Ai(iV = 2L) = 0 but A 2 {N = 2L) > 0, 
signaling a TI for f//f < Ucr- 

The MI extends for U/t S> 1, where a perturbative expansion shows that the system 
can be described by a spin model in a paramagnetic phase: under the assumptions 
U J,t and p = 2, one atom is trapped in each site of the two legs. We thus introduce 
the Pauli operators {i = x,y,z) acting on the local effective Hilbert space spanned 
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by the two spin states a = ±1 of the atom located at the site r of the chain r = ±1. 
We obtain the following second-order perturbative spin Hamiltonian: 

Hpert = ^ ^ '^eSVr,TVr+l,T “I" COS B (j]r^TVr+l,T “1“ Vr,TVr+l,T) 

r^T 

+ Jeff sin B - f]lr€+l,r)] + Kesvi,r=l€,r=-1 (10) 

r,j 

where Jeff oc t^/U and iJeff oc J^/U. In this regime, the term proportional to H 
dominates and the ground state of Jfpert is close to a trivial product state in which 
all the spins are oriented in the x direction. We expect that such state, characterizing 
the MI phase, might be adiabatically connected to the trivial band insulator at H = 0 
and J > Jc,i. Our numerics does not suggest the existence of a further phase transition 
between the Mott and the trivial band insulating phases. 

The phase diagram in Fig. shows that the topological region appears clearly as a 
thermodynamic region, within a well dehned parameter regime. We emphasize that, for 
1-8 ~ J ~ 2.2, the system is in a topologically trivial phase for H = 0, and enters the 
symmetry-protected topological phase when the interaction parameter U is increased. 
Therefore the interaction is not necessarily obnoxious to the purpose of experimentally 
obtaining the topological phase but, on the contrary, it can also drive the system into 
it by shifting the position of the critical point. This has been verihed also in the 
corresponding topological superconductor systems [351 EHl EH, where the addition of 
repulsive interactions is proven to expand the topological phase for certain ranges of the 
physical parameters (see also [3^ 140] for related models in terms of Majorana modes). 
This means that, for some particular value of J > Jc,i the presence of a repulsive 
interaction allows the formation of edge modes otherwise absent. A similar behavior is 
also observed in 2D systems with time-reversal invariance [66] • Let us stress, however, 
that this has nothing to do with the physics of fractional Chern insulators, where 
interactions drive the system into distinct (strongly-correlated) topological phases. As 
the phase diagram clearly shows, there is only one TI phase, which is strictly equivalent 
to that of the non-interacting system. Importantly, the phase diagram in Fig. [^ shows 
that interactions have a non-trivial role in tuning the system in and out the TI phase. 

In order to further clarify this last point, we now investigate in more detail the 
properties of the topological phase in the interacting system. Numerical investigations 
reported in Fig. show that the signatures of the non-interacting TI persist in the 
presence of interactions. First, the density prohle of the gas allows for a clear 
identihcation of the presence of fractionalized edge modes located at the boundaries 
of the ladder via the computation of {hf) and (h*). Indeed, for U/t < Ucr, Fig. 
shows that (h*) saturates to 0.5 within few sites, which is strongly different from the 
behaviour for U/t > Ucr- It is interesting to observe that within the topological region 
the localisation length of the edge modes has a weak dependence on U/t. Second, the 
system displays also within the interacting region a non-zero winding number associated 
with the second band of the system. As in an interacting system bands are not well 
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Figure 8. Properties of the interacting topological phase for J/t = 1, Vl/t = 1.8, 
B = f all the analyzed systems are on the line J/t = 1 in the phase diagram in 
Fig. i L = 72 and MPS bond dimension D = 200. (top, left) Density profiles {hf} and 
(np of the gas with U/t = 1 for N = 2L and N = 2L ± 1. The data clearly show the 
presence of localised and fractionalised edge modes, (top, right) The plot of \ {nj) — 2\ 
highlights the localisation of the edge modes, (bottom, left) Spin winding associated 
to the second band of the model S 2 for several values of the interaction, within and 
without the TI. (bottom, right) Entanglement spectrum (60 largest eigenvalues) for 
N = 2L and U/t = 1. 


defined, the winding number is computed by subtraction of the data relative to p = 1 
to those relative to p = 2 (see similar discussion in Sec. |^. This robustness of the 
spin winding number against local interaction is consistent with similar results in two- 
dimensional systems [H]. Finally, on a more abstract side, the analysis of the Schmidt 
spectrum presents the robust two-fold degeneracy of symmetry-protected topological 
phases |SI|. 

5. Physical realization of the model 

The physical realization of the ladder system in Eq. ([^ can be obtained by extending 
the 2D setup elaborated and realized in Ref. [1]. The present proposal builds on a 2D 
optical superlattice subjected to a well-designed time-modulation as displayed in Fig. 
Along the y direction a superlattice potential is used to partition the lattice into a ID 
array of isolated ladders. Hopping between the two legs of the ladder corresponds to 
transitions = —1 -H- -fl, see Fig.[^ The main challenge in realizing the Hamiltonian in 
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Figure 9. Schematic drawing of the proposed experimental setup. The 2D lattice 
configuration consists of a spin-independent superlattice along y to isolate individual 
ladders from each other and a superlattice potential along x, which creates a spin- 
dependent energy offset Atr^ between neighboring sites in order to inhibit tunneling. 
Tunneling is then restored resonantly with two pairs of beams denoted as r and b 
following the scheme introduced in Ref. [4j. Each of the pairs consists of a standing 
wave along x and a running-wave along y. For uj^/b = <^i,r/b — W 2 ,r/b = iA and 
Qr/b = ^i,r/b ~ ^2,r/b = (1; 1) ’ 7r/(2a) an effective flux d* = tt is realized with spin- 
dependent complex tunneling-matrix elements, here a is the lattice constant of the 
potential along x. 


Eq. ([^ consists in engineering the spin-dependent complex matrix elements for tnnneling 
processes taking place along the legs of the ladder. In the following, we will show that 
this can be achieved by combining a spin-dependent snperlattice potential x (Fig. [^, 
which introdnces a spin-dependent energy offset A between neighboring sites, together 
with the space-dependent time-modnlation of the lattice discnssed in Ref. [1]. 

We start by considering the time-independent part of the system, which can be 
described by the 2D tight-binding Hamiltonian 

Hq Jx ^ ^ ^^m+l,n^rn,n “h H.C.^ -|- (H) 

m^n 

~ '^y ^ ^ (^^m,n+l^rn,n + H.C.^ -|- — ^^( — 1) d\^ .^(Jzdxn,ni 

m,n m,n 

where m and n label the horizontal and vertical integer coordinates. The spin-dependent 
staggered potential conld be realized, for instance, by considering an appropriate anti¬ 
magic wavelength, for which the polarizability is opposite for the two spin species [5T] . 
In order to keep the bare tnnneling processes of strength ^ spin-independent the 
remaining lattice potentials need to be created nsing a magic wavelength, for which the 




















Methods for detecting charge fractionalization and winding numbers... 


17 


polarizability is the same for the two spin-species. 

The bare tunneling is suppressed along the legs due to the offset A which 

allows for a complete control over induced-tunneling-matrix elements, such as those 
realized by modulating the lattice resonantly [53]. Following Ref. [1], the modulation is 
taken to be produced by two pairs of laser beams with frequency difference iWr/b = ±A 
(Fig, g in order to restore resonant tunneling. The corresponding time-dependent 
potential dehned by these four lasers is then of the form 


Vi't) = 


^m,n^rn,n 


[n(m, -1-n*(m, n)e , 


( 12 ) 


with the resonance condition ca = A, and we choose the laser phases in such a way that 


v{m,n) = - |cos 


TT/d) ^ ^ gmn+ii?/2| 


/• 


(13) 


This requires a stabilization of the phase of the modulation relative to the static 
lattice potential, which is challenging and was not yet demonstrated in previous 
realizations HHUSni El Sj. This specihc choice of the potential v{m,n) is made in 
order to independently address successive hopping terms along the x direction, which is 
generally required when engineering Peierls phase-factors in superlattice structures, see 
Refs. [53l El and below. 

The time-evolution of the system is ruled by the Schrodinger equation idtf) = 
H{t)ip., where H[t] = Hq+V (t) is dehned by Eqs. 0 and ( [I^ . The long-time dynamics 
of the system can be suitably described by an effective-Hamiltonian approach |H3] , which 
is valid in the high-frequency regime u ^ oo. Since the static Hamiltonian Hq contains 
a staggered-potential term that explicitly diverges linearly with A = ca, we hrst apply 
the unitary transformation 


A 


= R{t)fj = exp i-iWtj f), W = — 


(14) 


which removes the diverging term. The effective Hamiltonian can then be derived in 
this moving frame, using the method of Refs. [531 EH] (see also Ref. [55]h 

For the sake of simplicity, let us hrst consider the dynamics associated with the 
species = -|-1. For these atoms, the transformed Hamiltonian reads: 


H(t) = R^(t) ffo + V(t) R(t) - W = + V-e 




(15) 


where 


m^n m odd,?! 

^ ^ y ^ "hm^nV {rn, n) — Jx y ^ • 


m,n 


m even,n 


( 16 ) 
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We describe the time-evolution of the system dictated by H{t) through the evolution 
operator, which we partition as 




U{f) = 


(17) 


where the effective Hamiltonian ffes describes the long-time dynamics, and where the 
operator K{t) captures the micro-motion. Following Ref. [63], we find that the effective 


Hamiltonian associated with the general time-dependent Hamiltonian in Eq. (15) is 
given by: 


(jJ 


(18) 




u 


®L+l,n«m,n (m -|- 1, u) — v{m, n)) + 


m even, n 


-f- {v{m -l,n) - v{m, n)) 


+ H.C. 


JxK 

\/2u 

JtK 


E 


aL^i cos 


‘'m+l,n'^rn,n ' 


-l,n^rn,n COS ^ 6 




m even, n 


+ H.C. 


V2i 


UJ 




I2 


E 

m even 

E 

m odd 


®2m-|-l,n®2m,n “1“ ®2m,n®2m—l,r 


'^2m+l,rfl2m,n ®2m,n®2m—l,n 


+ H.C. 


The irrelevant sign change in the tunneling matrix elements (i.e. in the last line of 


Eq. (18)), can be removed by applying an additional gauge transformation 


G = exp 


'y ^ ®2m,n®2m,n 


m odd,n 


= G\ G" = 1, 


(19) 


which indeed reverses the sign of the tunneling terms a 2 m+in® 2 m,n and a 2 mn® 2 m-i,n for 
m odd only. In this way, the final effective Hamiltonian describing the dynamics of the 
az = +1 species yields 


HeS —t GHe^G — 

^ m,n 


-t A m7r+iB/2 i pr „ 


( 20 ) 


which is indeed the tunneling term in Eq. (A.l) for = -|-1 atoms. In the case of the 


az = —I species, the staggered potential is reversed, so that the even and odd sites must 
be inverted. This results in the final effective Hamiltonian: 


= V E + H.c.) , 

m,n 


( 21 ) 
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where (—1)"' is equivalent to the operator in Eq. ( A.l[ ). 

The Zeeman term present in the Hamiltonian Q can be directly generated 
by two resonant coupling potentials V/^^(t) = Xlm n cos(r'it)am„nO'a;am,n- Indeed, 
considering the bare atomic frequency 3> A between the two sublevels, they are 
effectively separated by the position dependent energy offset + (—1)'"A created by 
the spin-dependent potential in Eq. Q' Under the transformation R{t), the effect of 
the coupling becomes 


R\t)Vi^^{t)R{t) = ^ 2Qe cos (uit) d\^^^a+dm,n + H.c. 

m,n 


( 22 ) 


where ct+ = {ax + iay)/2. This terms commute with the gauge-transformation operator 
G and, by choosing the frequencies 012 = ± A, we recover the required Zeeman 

term ilax in Eq. Q, both for even and odd sites, through the standard rotating-wave 
approximation. 

We note that the effective Hamiltonian Q was derived at first order in 0 ;“^, in a 
basis provided by two commuting unitary operators: R{t) and G. It is important to 
notice that these latter operators commute with the r operators, so that they neither 
affect the static hopping term Jtx, nor the static potential difference /ir^: this indicates 
that these static terms can be directly included in the (effective) Hamiltonian ([^. The 
latter remark is also valid for the spin-independent potential /xq- Therefore, we conclude 
that the application of two pairs of Raman lasers and a radio-frequency field, combined 
with the spin-dependent staggered potential directed along the ladder, allows one to 
generate all the spin-dependent terms in the ladder Hamiltonian ([^. 

Importantly, we emphasize that the Hubbard interactions are also unaffected by 
the aforementioned transformations G and R{t). Thus, the effects of interactions can 
be directly incorporated into the effective Hamiltonian ([^, at first order in uj~^. 

Finally, the time-evolution operator in Eq. (0 is then completely determined by 
computing the kick operator K{t), which is readily calculated using the expression 

[SSI ES] 


k{t) = ^ 


lU 


U+ 


Au)t 




V-e 


2k, 


~ — y^hm,n\v{m,n)\sm{ujt + 9m,n)y (23) 

mfa 


where we assumed that k^ Jx and we defined 6m,n = arg[u(m,n)]. 


6. Conclusions 

In this work we presented a ladder setup for ultracold fermions subject to both the 
presence of an artihcial tt-Aux magnetic potential and a spin-orbit coupling. Such a 
model may be seen as the synthesis of two accessible experimental techniques to realize 
synthetic gauge fields in optical lattices: on one side, the realization of complex tunneling 
matrix elements using time-modulated optical lattice [3lHJISI091 [50] , and, on the other, 
the implementation of spin-orbit terms through the control over internal atomic degrees 
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of freedom EE]. Analogously to models already discussed in the context of nanowires 
[29l |30] , the combination of these two elements gives rise to a particle-hole symmetry 
which protects non-trivial topological phases within the Bogoliubov-de Gennes and 
chiral classes of topological insulators and superconductors. As the model conserves the 
number of particles and can be realized with state-of-the-art experimental techniques, 
our results might give a substantial advance towards the observation of topologically- 
protected zero-energy modes in fermionic systems with and without interactions. In 
particular the physical realization that we present does not require the engineering of 
any pairing mechanism, neither the coupling to external molecular gases or superfluids, 
as, for example in [HI [HI IHl IH], nor an interchain pair-hopping, as exploited in the 
ladder model presented in na. 

No interaction is indeed necessary for the appearance of the symmetry-protected 
topological phase. Consequently, the realization of the particle-hole symmetry through 
the ladder geometry is not as robust as its counterpart in topological superconductors 
due to the absence of a true superconducting gap, thus it should be considered an 
extrinsic feature. However, due to the absence of disorder and to the high degree of 
isolation in ultracold atom setups, we expect the topological features of the system to 
be experimentally detectable. For example, the introduction of a trapping potential does 
not spoil the observation of edge modes, despite breaking the particle-hole symmetry. 

We analyzed two experimentally relevant signatures of the appearance of topological 
phases: the presence of fractionalized edge modes, detectable through site-resolved 
density measurements (as recently reported in Refs. [59l EOl El]), and the winding 
behavior of the spin degree of freedom, which can be observed through spin-resolved 
time-of-flight imaging. In particular, we have shown how to detect these observables even 
in the presence of a trapping potential, which induces soft boundaries, often believed to 
be particularly disruptive for the detection of topological signatures. 

Our study has also considered the effect of a contact repulsive interaction; this 
is possible thanks to the presence of the spin degree of freedom that differentiate the 
main features of our model from its spinless counterparts as the SSH and its interacting 
generalizations (see for example mm)- Apart from mapping out the phase-diagram 
of the model, which entails two gapped phases, with and without topological properties, 
we have found that a Hubbard interaction can enhance the extension of the topological 
phase, instead of being detrimental. Furthermore, the spin winding number introduced 
in the article provides a good topological parameter also in the interacting case, where 
the usual order parameters based on single-particle wavefunctions fail. 

Concluding, we mention that the model under scrutiny may be an interesting 
platform for the study of further fractionalization effects, based on a particular fine 
tuning of the parameters, which can be reminiscent of the physics of parafermionic 
zero-energy modes (see Ref. [69] for a recent review) in the spirit of Ref. PI. 


Methods for detecting charge fractionalization and winding numbers... 

Acknowledgements 


21 


We acknowledge enlightening discnssions with Rosario Fazio, Fabian Hassler, Sylvain 
Nascimbene, Lnca Taddia and Andrea Trombettoni and thank Simone Barbarino for 
carefnl reading of the mannscript. We thank Davide Rossini for providing the nnmerical 
algorithms based on Matrix-Prodnct States employed in this work. L. M. acknowledges 
snpport from the Italian MIUR throngh FIRB project RBFR12NLNA and Regione 
Toscana FOR FSE 2007-2013. M. A., H.-H. T. and M. B. acknowledge snpport from 
the EU grant SIQS. N.G. is hnanced by the FRS-FNRS Belginm and by the BSPO 
nnder the PAI project P7/18 DYGEST. 


Appendix A. Symmetries of the model, order parameter and spin winding 
number 


In this Appendix we examine the non-interacting ladder model and we discnss, in 
particnlar, the relation between the order parameter W, that distingnishes trivial and 
topological phases, and the observed spin winding nnmber. 

The Hamiltonian Q is translationally invariant and can be also expressed in 
momentnm space as 77 = d\T-i{k)dk where: 


77(7) = 2tTz cos(7?/2) cos k — 2tTz(Jz sin(R/2) sin k -|- + Jt^ + -|- / iq . (A.l) 

We observe that its kinetic term corresponds to a Peierls snbstitntion k ^ k + 7^. 

Its spectrnm is symmetric for a tranformation mapping B ^ B + 2ti and t —)■ —t. 
Therefore, even if the Hamiltonian is periodic in B with period 47r, we can restrict 
onr stndy to the case 0 < 7? < 27r. Besides, we can consider only positive valnes of 
G and J becanse their sign trivially depends on the chosen basis for the spin and 
psendospin. In particnlar 77(7,—G) = azT-L{k,VL)az and 'H{k,—J) = TzB^k, J)tz. 
Analogonsly we consider only 7 > 0, becanse of the symmetry between the two chains 
77(7, —t, —/i) = Tx'H{k, t, ij)tx. Hereafter we rescale all the energies in nnits of t in snch a 
way that, below, we will always consider 7 = 1. Finally we observe that the expectation 


valne of ay is always nnll becanse the Hamiltonian (A.l) is real. 


The system is characterized by an anti-nnitary time-reversal-like symmetry T = a^ 


T'H{k)T^ = W{-k), 


(A.2) 


and, for /tq = 0 (the system is exactly at half hlling), we obtain the particle-hole 
symmetry C = Tyay\ 

Cn{k)C^ = -n*i-k). (A.3) 

These non-nnitary symmetries characterize the topological symmetry class BDI (see, 
for example, [211 123 ]), which is also characterized by the nnitary chiral symmetry 
P = TC = Tyaz- 

Additional terms in the Hamiltonian may break the C symmetry which is indeed 
fragile, dne to the lack of a physical pairing interaction; we emphasize however that 
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such perturbations become significant only if their magnitude is comparable with the 
energy gap. To this concern, in a cold atom gas, the presence of noise and defects are 
negligible and the main effect we must consider is the trapping potential bringing to a 
space dependent chemical potential /iq. As discussed in Sec. however, the local shift 
in energy provided by the trapping allows in general to isolate topological regions of 
the chain with the effect of binding fractionalized modes at the interface between these 
regions and the trivial ones. 

An additional Zeeman term proportional to cXx, which breaks only the time 
reversal symmetry, brings instead the system in the symmetry class D, which is still 
topologically non-trivial in one dimensions. Therefore, this sort of term does not alter 
in a fundamental way the properties of the system. 

One-dimensional systems in the BDI class possess topological phases labelled by a 
topological invariant in Z [23]. This topological invariant, can be evaluated by exploiting 
the chiral basis defined by the symmetry P [56|. In this basis the Hamiltonian assumes 
the simple form: 

( 0 A.{k)\ A A/i + J + 2icos (ffc) iVL 

(A.4) 

If det(A(/c)) 7 ^ 0 for all k then the system is gapped. In this case, detA(fc) = 
|det A(fc)| and the winding number of f{k) constitutes the topological invariant 
which distinguish topological and non-topological phases [SS]- In this model, this 
winding number may assume only the values 0 or ±1 characterizing trivial and 
topological phases respectively. In more detail, one obtains: 



det A{k) = — — /i^ -I- — 2 cos{B) — 2 cos(2A;)-|- 

— 4/icos(H/2) cos(/c) -I- 4iJsin(i?/2) sin(fc). (A.5) 


This determinant is purely real for /c = 0, tt or for B = 0. For B = 0, its phase f{k) 
cannot change its value from 0 to tt unless det A crosses zero and the gap closes; thus 
the case H = 0 is either trivial or gapless. 

In all the other cases, for 0 < H < 2n, the topological invariant W is evaluated by 
considering the behavior of f between k = 0 and k = n. In particular, the tangent of 
the phase f is given by: 


tan [^(/c)] 


4 J sin ^ sin k 

— 4/i cos A cos k — 2 cos B — — 2 cos(2fc) — /i^ 


(A.6) 


where the numerator is always positive for B G (0,27r) and k G (0,7r), and always 
negative for B G (0,27r) and k G (tt, 27r). Therefore, going from /c = 0 to /c = tt, ^{k) 
must be always included in (0, vr) since sin(^) > 0 in this regime, whereas for k which 
goes from vr to 27r, f{k) must be either in (—vr, 0) (if f{k = 7i) =0), such that its winding 
number vanishes for k —)■ 27r, or ^{k) G (vr, 27r) which implies a final winding number 
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Figure Al. Different behavior of the phase ^ as a function of the momentum k for 
a trivial phase (red dashed line) and the topological phase (blue line). In the trivial 
case the values of ^ at A: = 0, tt, 2tt are equal and the winding number parity is W = 1; 
for the topological phase, instead, = 0) and ^{k = 27r) differ by 27r. The phases 
were calculated for /i = J = 1 and B = 47r/3 for D = 1.2 in the topological region and 
D = 2.3 in one of the trivial regions. The shaded regions are forbidden for B G (0, 27r). 


equal to one (see Fig. Al). Thus the parity of the winding number results: 


W = sign 


— 4/tcos(i?/2) — 2 cos(i?) 


(A.7) 


Q? — 2 — — jj? + Ajji cos(i?/2) — 2 cos(i?) 

W = 1 when the winding number is 0, whereas W = — 1 for the winding number 
being ±1. Therefore, in terms of the parameter 12, two phase transitions appear at: 


and 


12^ ^ + min 


12 ^ 2 = '^^ + 


/i 


B 

2 cos — 
2 




B 

2 cos — 
2 


B 

fi — 2 cos — 


B 

B + 2 cos — 


(A.8) 


(A.9) 


The system is in a topological phase (W = —1) for 12^^ < 12^ < 12^2 whereas for 
12^ < 12c;i and 12^ > 12^]^ we obtain trivial phases (W = 1). At the transition points 
12^ = 12 ^^, 12^25 fhe gap closes respectively for fc = tt, 0, consistently with the results 
[301 [70]. The corresponding critical values for J are given by: 


m 


Ji 2 - =F 2cos(i?/2))' 


(A. 10) 


Depending on 12,/i and B there can be 0,1 or 2 phase transitions as a function of J. 

In order to evaluate the spin expectation value in the plane h — z, it is useful to 


adopt the Hamiltonian form in Eq. (A.4) which allows a simple diagonalization. In this 


basis the components of the physical spin, and become: 


= a. 


r, 


y 1 


S'" = Tq. 


(A.ll) 
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where cTj and r* are Pauli matrices in this basis. To obtain the eigenstates of the 
Hamiltonian we consider: 




A{k)A\k) 0 A 
0 A^{k)A{k) J 


(A.12) 


where, AA'^ = (A'^A)* since A 
this basis is: 


A^. Therefore the generic form of the eigenstates in 



(A.13) 


where if are the two eigenvectors of AA^ . In particular, we are interested in the 
eigenvector corresponding to the lowest eigenvalue of AA\ since it determines the second 
and the third bands of H which are topologically non-triviah The value of a is fixed by 
the equation: 

Alp* = ±£e"*“ (A. 14) 


where ±£ are the energies of the two intermediate bands of the system. 

The expectation value of in the second band does not depend on the phase a 
and it is easily evaluated. It results: 


sAk) ^ {s^) 


2 sm{k) (sin(i?) cos(fc) + /isin (^)) 
N{k) 


(A.15) 


where N{k) is a positive quantity. In particular we must distinguish two cases: 
|/i| ^ |cosi?/2|. 

For fi > 2|cosH/2|, S'! > 0 for A; G (0,7r), thus S'! behaves like the imaginary 
component of the phase exp[i^(A:)] and it can be showed that the sign of S'! at /c = 0 
and k = n is equal to the sign of the denominator of Eq. (A.6), therefore the spin 
winding number of the second band and the winding number of the phase fA) coincide. 

For 0 < /i < 2 |cosi?/2|, S^k) has a further zero between k = 0 and k = n, therefore 
the spin vector is aligned along x in three different points in the interval k G [0, tt]. This 
can be seen in the last column of Fig|^ (obtained for /r = 0.5 and B = 7r/2), where the 
blue line, representing the average value of the spin in a system without the trapping, 
always crosses the horizontal axis three times. The spin winding number S 2 is null if 
and only if S^ik) has the same sign in these three points (first and last rows). However 
if S^ik) assumes the same sign for /c = 0 and fc = tt, then the total spin winding number 
for k going from 0 to 2t: is even, due to the time-reversal symmetry T. Otherwise it is 
odd. At /c = 0 and k = tt the sign of Sipik) assumes, respectively, the same value of the 
signs of the numerator and denominator of W in Eq. (A.7). Therefore the parity of the 
spin winding number, (—always coincides with W. 

Concerning the behavior at the intermediate zero of S'!, 


kc = arccos (—/isin(i?/2)/sin(i?)), 


(A.16) 
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one has S^ikc) = cos [a{kc)] oc fl — J, therefore in the trivial region where 2 the 

spin winding nnmber is trivial. The other trivial region, is instead divided 

into two snbregimes characterized by a spin winding number ^2 = 0 or 1 S 2 = 2 separated 
in J = n. This is shown in the last column of Fig. where the blue lines either does 
not wind around the origin (first row, corresponding to hi < J < 12 ^, 1 ) or winds twice 
around it (second row, hlcp > f2 > J). 

Similar conclusions hold for negative values of /i and it can be shown that the 
behavior of the spin winding number is not affected by open or periodic boundary 
conditions. 
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